Peculiarities of anisotropy and polarization as an indicator of noises in the CMB maps 
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We discuss some new problems of the modern cosmology which arose after the BOOMERANG 
and MAXIMA- 1 successful missions. Statistics of high peaks of the CMB anisotropy is analyzed 
and we discuss possible inner structure of such peaks in the observational data of future MAP 
and PLANCK missions. We have investigated geometrical and statistical properties of the CMB 
polarization around such high isolated peaks of anisotropy in the presence of a polarized pixel noise 
and point sources. The structure of polarization fields in the vicinity of singular points with zero 
polarization is very sensitive to the level of pixel noises and point sources in the CMB maps. 
PACS number(s): 98.80.Cq, 95.35.+d, 97.60.Lf, 98.70.Vc. 

I. INTRODUCTION 

Observational data by BOOMERANG and MAXIMA- 1 [1,2] open a new epoch in the investigation of the CMB 
power spectrum at the large multipole numbers I. The measured angular power spectrum shows a clear peak at 
angular scales corresponding to the spherical harmonic multipole number I « 200. This is the Sakharov's peak. 
However the structure of the CMB anisotropy power spectrum at I > 400 is still unclear and next generation of 
satellite experiments (such as MAP and PLANCK) are needed. 

There exist a few factors which play a leading role in the future experiments. Namely, both the MAP and PLANCK 
missions will provide significantly greater percentage of sky coverage than BOOMERANG and MAXIMA-1. In 
addition in the PLANCK mission there are two HFI channels ~ 545GHz and i' ~ 857GIIz which provide FWHM — 
5arcmin resolution. Future polarization measurements are also very important. 

In contrast to BOOMERANG and MAXIMA-1 in both high frequency channels of PLANCK contamination fore- 
grounds are of a prime consideration [3]. Of course, there is a wide variety of "technical" methods which could be 
applied to obtain a well-cleaned data over a wide range of angular scales. For example, one can use the frequency 
dependence of the various foregrounds (dust emission, synchrotron, thermal emission, point sources etc.) [3,4]. In 
practice this method meets serious difficulties because of the angular distribution on the sky and the frequency 
dependence of the foreground components are not well known. 



That is why we would Hke to focus our attention on the investigation of local characteristics of the CMB anisotropy 
and polarization. 

In this paper we analyze the statistics of the peaks (maxima and minima) in the CMB anisotropy maps. We compare 
this statistics for the BOOMERANG and MAXIMA-1 maps with the statistics for the future MAP and PLANCK 
observations and predict some properties of the peaks and their shapes in the future PLANCK observations. The 
main point of our attention is the analysis of the structure of the polarization filed around singular points with zero 
polarization. We propose to use the results of the analysis for the estimation of the level of noise and foregrounds in 
the maps. 



II. STATISTICS OF PEAKS IN THE CMB ANISOTROPY MAPS 



In this section we compare statistics of peaks (maxima and minima) for the BOOMERANG and MAXIMA-1 maps 
with the statistics for the future MAP and PLANCK observations. 

Let us consider the model which is close to the real situation provided by the successful MAXIMA-1 balloon mission. 
According to [2] in the MAXIMA-1 pixeled map one can find a high amplitude peak in AT/T distribution on the 
observed region of the sky with following coordinates: declination ~ 58.6 degrees; RA= 15.35 hours. In the Wiener- 
filtered map this peak has an amplitude AT ~ 2.3 2.5a and it decreases monotonically down to la cross-level at 
15.2 < h < 15.4 hours, 58.5 < 6 < 60 degree. 

Below we predict some definite properties of the structure of such peaks in the future PLANCK CMB measurements 
performed with higher angular resolution than the current measurements. Analogous prediction has been done in [5] 
for the Tcncrifc Experiment using COBE DMR data although with the different technique. For example, we show 
that more accurate measurements will not reveal the inner structure of the peak in a form of new high (> 1.5a) 
peaks inside the area mentioned above. We assume that the fluctuations of the CMB anisotropy (or polarization) are 
realization of a random Gaussian process. It is well known that for such random processes all statistical properties of 
a signal are determined by the correlation function of the anisotropy 

Crie) = ^ 5^(2^ + l)Ci,TW{l)Pi{cose), (1) 
I 

where Ci^t is the power spectrum of fluctuations of anisotropy, W{1) is the window-function which depends on the 
strategy of the experiment, Pi{cos9) are Legendre polynomials. Analogous expressions can be written for Q and U 
components of polarization: Cq and Cu- 

To investigate topology of the anisotropy and polarization maps we deflne the spectral parameters (as in [6,7]) 

= (-ir*!2^^^; c. = 2sin^; (2) 



^V2^; 

a2 



9cf = -C{0)/ 



where « = 0, 1, 2, and 9a is the antenna beam. 

In the real experiments the antenna has a finite resolution and the spectral parameters depend on the antenna 
beam 0a and therefore on the number of Doppler peaks which could be resolved by the antenna. This means that in 
the future PLANCK maps the structure of the high AT/T peaks can be different from the corresponding structure 
in the BOOMERANG and MAXIMA-1 maps. 

For all cosmological models the power spectrum of the CMB anisotropy Ci can be described in terms of a sum of 
the Gaussian peaks centered at the points of the maxima Z„ , (/ > 30) 
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where n is the number of the peak, d„ is the width of the peak, is its position, is its amplitude and N is the 
normalization factor for low multipole range (for example, the COBE data normalization). The last term in Eq.(P) 



^We will omit the index T in the subsequent discussions 
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accounts for the Silk-damping at the angular scale s. Note that we do not include in Eq.(|^) any low and high multipole 
filters (the window function and beam) . These means that Eq.(l) describes the initial power spectrum of fluctuations 
on the sky without any smoothing. In reality both last factors are extremely important and their influence on the 
AT/T maps plays a crucial role. The power spectrum in the form of Eq.(^) after substitution into Eq.(||) gives 
information about the influence of each peak on the topology of AT/T maps. Thus, using the approximate Eq.(|^) 
we can investigate the influence of the first, second and subsequent Doppler peaks on the spectral parameters of the 
future maps of the MAP, PLANCK and other missions. 

Next point is connected with the window function W{1) of the experiments. In a small angular approximation we 
will model general properties of W{1) as follows: 
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where m — 2 for the single difference and m — 3 for the double difference scheme of the low multipole filtration. The 
exponent in Eq.(^ describes the antenna beam with 9a — 7.45 x 10~^ For description of the asymptotes 

of the multiplier in Eq. (0) we can introduce a function which matches both limits 
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where R ^ 0.03 is the characteristic angular scale at the low multipole filtration [see Eq. 
spectral parameters in Eq.(0) are: 



)]. So, in such a model the 
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Note that for the second and higher Doppler peaks we have I'^/d^ ^ 1 and only for the first Doppler peak l\/d\ ~ 5. 
For the analytical approximation of the integral in Eq.(^ we require the asymptotic If^/d"^ 3> 1 for all peaks in the 
power spectrum Eq.(0). Using this approximation we obtain the following result for the spectral parameters af: 
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where C is the Euler constant, = 6*^ + s'^, $(a;) = 2/y^ dxe ^ is the probability integral. As we can see from 
Eq. only the first Doppler peak is important for calculation of the variance ctq . The influence of the second and 
higher peaks is practically negligible due to the decrease of amplitudes An and d„/Z„. However these peaks determine 
(see Eqs.(10)-(12)) the topological structure of AT/T maps (for example, the number of maxima and minima at 
different thresholds f„cro — AT/T). Using Eqs.(8) and (9) we describe the realistic model at ^2^2 <^ i ^nd ^2^2 < i_ 
In such a model the density of all peaks at 1/ g (—00, 00) has an especially simple form 
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where Np^ and Np^^ are the densities of all maxima and minima correspondingly. The density of saddle points (of 
arbitrary hight) is 



Nsad{.-oo) = 2N+j^ (stcr) 



(11) 
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Let us now discuss a model in which all Doppler peaks are smoothed {A„ 
parameters 6** and 7 are the following 



= 0). In such a model the spectral 
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and the densities of all maxima and minima (with arbitrary hight) are Np^ = Npj^ = ^^^ i ^- If the sky coverage 
for some AT/T - experiment is fsky (for example, fsky — 0.3% for the MAXIMA-1 experiment), then the number of 
maxima (or minima) in the observational map is 



N„ 



16 



fsky 

0.003 



1° 



(13) 



According to [1,2] for MAXIMA-1 and BOOMERANG the antenna beam corresponds to FWHM ~ 10'. This 
means that in the absence of the Doppler peaks in the power spectrum we can find 576 maxima on the corresponding 
maps. However one can find from Eq.(7)-(9) that the presence of Doppler peaks in the initial power spectrum changes 
the number of peaks in the map down to 271 for the mentioned above observational data. So, after our analysis we 
can conclude that the influence of the Doppler peaks leads to decrease of the number of hot and cold spots on the 
map by a factor ~ 2. This result consists with the data of the MAXIMA-1 and BOOMERANG maps. 

The next problem which wc would like to discuss below is the following. How sensitive is the topology of AT/T 
map to the amplitudes of the second A2 and third A3 Doppler peaks assuming that the amplitude and the position 
of the first one are known? To give answer to this question we compare 7(^2, A3) and Npk{A2, A3) for the following 
models. In the first model wc take the amplitude Ai corresponds to the data [1,2] and positions and widths of the 
next peaks in power spectrum are the following: li = 210 with the width di = 95, I2 = 580 with the width ^2 = 110 
and I3 = 950 with the width ds = 130. Corresponding plots are in Fig.l. 





FIG. 1. The dependence of the 7(0:, y) (left) and Npk{x,y) (right) on the parameters x = lO^A^/Ai (horizontal a:xes), 
y = lO^As/Ai (vertical axes). The numbers in the lines correspond to the values of 7(0;, y) and Npk{x,y). Note that x = 26 
and y = 46 corresponds to the MAXIMA-1 and BOOMERANG amplitude of the first, and two next Doppler peaks. 



The second model (see plots in Fig. 2) corresponds to a hypothetical situation when the amplitude of the first Doppler 
peak is two times smaller than in the previous case. As one can see from Eq(7)-(9) and Fig.l, Fig. 2, in this second 
"toy" model the structure of the spectral parameters 7(^2,^3) and Npk{A2,A3) changes drastically. The number of 
the maxima increases up to more than 420 while 7 parameter conserves practically the same value: 7 ~ 0.4—0.47. This 
result is important for analysis of the global and local topology of the maps. In the BOOMERANG and MAXIMA- 
1 experiments the amplitude and position of the first Doppler peak in C;-power spectrum are measured with 10% 
accuracy. That means that theoretical predictions of the number of the CMB peaks in the observational maps could 
change from 263 to 279 due to this 10% imccrtainty. The difference in 16 peaks corresponds to ^ ^ N~^^'^ statistical 
fluctuations of the number of peaks A'^ in the map practically without any essential changes of the 7-parameter. For 
other aspects of the distribution of the peaks in the maps see [7-9] . 
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FIG. 2. The same as in Fig.l but for a toy model with the amplitude of the first Doppler peak is two times less than in Fig.l 



In this section we consider how the structure and shape of the high peaks in the observational maps change when 
one goes to higher angular resolution maps. In this case one can measure the structure of the AT- field around a high 
peak with more details. 

Can the future measurements reveal internal structure of the peaks which are found by BOOMERANG and 
MAXIMA- 1?. For example can the future measurements reveal new peaks in the fine structure inside the region 
around the ~ 2 - peak down to the v = 07 And if yes, what would be a typical height of such peaks? The answers to 
these questions depend on the peak-peak correlation in a high-resolution map. Two important properties of a random 
Gaussian field play a leading role in this problem. First, in the vicinity of a high maximum of a Gaussian field the 
random character of AT{x,y) distribution is broken, and the shape of the peak is regular. The typical scale of such 
a regularity is, of course, 9^ (see Eq.(^)). Second, as we mentioned above, the fractional variance in the peak number 
from one realization to another over a grid of area Qp is related to the peak-peak correlation function Cpk-pk [6,9]: 



where (N^f^) — n'^i^{i/t)^p] ^pki^t) is the integrated number density of the maxima with height v above some threshold 
vt- Note that the first term in Eq.([l^) corresponds to the Poissonian distribution of peaks. Recently Heavens and 
Sheth [9] performed analytical and numerical calculations of the peak-peak correlation function and showed that 
Cpk-pk goes to zero at 9 < 9^,, and reaches negative value Cpk-pk = — 1, at = 0. This result reflects the fact that 
the different high peaks cannot be located close to each other and, for example, two high peaks with the amplitudes 
i/j^ ~ ~ 2 2.5(7 should be separated by the distance 9 ^ 9^,. According to [9] the typical angular scale 6** for 
the most preferable ACDM cosmological models is close to 20 arcmin. This scale is twice greater than FWHM in the 
BOOMERANG and MAXIMA-1 experiments and 4 times greater than in the PLANCK-mission. However, it is worth 
noting, that inside the mentioned above region around the high peak there can exist up to ^ 10 low amplitude peaks 
{v < 1) in a high resolution map. These peaks are extremely important for the proof of the Gaussian statistics of a 
signal in such a map. Thus we can conclude that the isolated 2 ~ 2.5a peaks which were found in the low-resolution 
BOOMERANG and MAXIMA-1 maps will reveal themselves as isolated peaks in the PLANCK map. 

Let us come back to the discussion of the high peak at the location 6 — 58.6 degrees; RA= 15.35 hours of the 
MAXIMA-1 map. The position of this peak practically does not depend on the higher angular resolution of the future 
PLANCK mission and its amplitude can be described as follows. Let us imagine that an ideal experiment with a 
(5-function antenna beam finds the highest peak in the AT map at the position {6; RA). An amplitude of such a peak 
measured in units of variance is: 



where cro(m) corresponds to Eq.(g) at 9a = and £, = s. We suppose for simplicity that the distribution of AT{x,y) 
around the point of maximum is Gaussian (more comprehensive consideration can be found in [6]) with the charac- 
teristic scales a and b = na, k is & constant: 



III. STRUCTURE OF THE HIGH PEAKS OF THE CMB ANISOTROPY. 
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where the parameter a is proportional to the typical correlation scale of the initial signal. Following [6] we can describe 
a local shape of the peak of height v, measuring radial curvature T and "ellipticity" e using polar coordinates and ^: 
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Let us denote 6 = + and cos 21^ = x^+y^ ■ 



After that we find: 
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In terms of eccentricity e = 2 [e/(l + 2e)]^^^ the last relation in Eg. (|l8|) can be transformed to 
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Let us consider two experiments which measure the same part of the sky around a peak with different angular 
resolutions 9i and 6*2. Suppose that 9i corresponds to the MAXIMA- 1 FWHM and 6*2 corresponds to the PLANCK 
(^1 ~ 2^2)- We denote an amplitude of the maximum in the low-resolution experiment by Vmaxima and for the 
high- resolution experiment by Vpianck ■ For such models an amplitude of the peak is 
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where index j — 1,2 corresponds to 9i and ^2, T{x,y) and r'{x',y') are vectors in the Cartesian coordinate system 
centered at the point of the maximum. The shape of ATj{x,y) from Eg. (p^) for the low- and high-resolution 
experiments is 
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This curve determines parameter ^2 = (52 ^ Q2y^^2 _|_ g2^ which we can measure in the vicinity of the peak at some 
threshold iyt'^Q^\ where ctq^'' is the variance of fluctuations in the low-resolution experiment. As a result the amplitude 
of the peak from Eg. (15) is 
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For the high-resolution experiment we obtain from Eg. (15) 
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and while the difference between a^™\ CTq^-' and CTq^'' is logarithmic only (see Eg.(^) we have 
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where /i = — < 1. 

^ a 

For example, for the MAXIMA-1 peak at the coordinates 5 = 58.6 degrees, RA~ 15.35 this ratio is ~ 1.2 1.4. 
Taking into account this result vi ~ 2 3 peaks in the MAXIMA-1 map will be transformed to higher peaks in the 
PLANCK map {1^2 — 3 ~ A), 9i ~ 26*2. This conclusions are important for our discussions in the next section. 
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IV. STRUCTURE OF THE POLARIZATION FIELD AROUND THE POINTS WHERE POLARIZATION 

OF THE CMB SIGNAL VANISHES. 



The structure of the polarization field and its auto- and cross-correlation with the anisotropy distribution on the 
sky has been studied in many papers (see for review [10]). Below we consider the structure of the polarization field 
around such singular points where P = in the CMB signal taking into account possible foregrounds and noise. This 
analysis can help in the estimation of a possible level of the foregrounds and noise in the future MAP and PLANCK 
polarization maps. 

Of course, measurements of the CMB polarization around the points where it vanishes are a very difficult observa- 
tional problem. However it is not necessary to go exactly to the point where + U'^ = (where Q and U are Stocks 
parameters), because the structure of the polarization field can be determined by the pattern created by the flux lines 
around peculiar points. We note also that in spite of the fact that these points are distributed almost randomly over a 
map it is more probably to find them in the vicinities of high peaks of the AT/T anisotropy. Indeed, in coming paper 
we will show that there is non-zero correlation between high maxima of anisotropy and singular points of polarization 
field. In this paper we consider the problem in the continues (non discrete) measurements of the Q and U values to 
reveal the main qualitative behaviors of the polarization field. The correspondingly consideration for pixelized maps 
see in [11]. We shall assume that in the observational map Q and U components of Stokes vector are the sum of the 
initial CMB signal {Q,U) and the Gaussian noise {nQ,njj) which model the influence of pixel noise and polarized 
point sources on the primordial CMB signal. We shall assume also that statistical properties of the noises uq and 
nu are identical and their correlation functions are Cf^ = Cf^ = (J^Sij, where is the variance. This means that we 

neglected the correlation between noise of different pixels. Under the assumption mentioned above we can expand Q 
and U near the points where Q = U = a,s follows: 

QKiqiX + (722/ /25^ 
U w uix + U2y, 

where x and y are Cartesian coordinates of the system centered at the zero-point of polarization; qi, q2 and ui, U2 
are the corresponding partial derivatives. Thus the modulus of the polarization vector = + U"^ is 

P^^P^ + 2(QnQ + Unu) + «q + n^, (26) 

where Q and U correspond to Eq.(|25|). 

Lets us firstly describe the model when uq and nu are negligible with respect to the regular functions Q and U . This 
situation takes place when + u'^ ^ Uq + nfj or P^ ^ a^, where cr^ is variance of the noise and in Eq. (^6|) we can 
neglect all terms excluding P^ . In such a case we obtain from Eq.(p6|): 
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ail = qfui, 022 ~ li^i'j 0.12 = 0,21 — qiq2 + U1U2. Rotating the coordinate system through the angle 

2ai2 
ail — a22 



tg2^ = — ^'^"„„ one can transform Eq.(p7|) to the "standard" form: 



P^ ^ Ax"^ + By'^ , (28) 

where A — ^ [an + 022 + Vl + t^(aii - 022)] ; B = ^ [an + 022 — + t^(aii — 022)] : A - B > 0; t —tg2ip and x and 

y are the coordinates after rotation. 

As one can see from Eq.(G8h the condition P — const determines the ellipses with a peculiar isolated point inside 

2 2 2 

them P" = 0. Therefore the structure of the field P^{x, y) around points P = for the case P ^ Uq -f is rather 

2 

simple. But when we come close enough to the point P = this condition will be violated because the level of noise 

is approximately constant > near the point P = 0. Now in Eq.(^6|) we take into account the second term but 
neglect ng +n'^, so 

P^ -P^ =2(QnQ + Unu). (29) 

As one can see from Eq. (|9|) the function f — P^ — P"^ is a nonuniform random function due to a regular character 
of Q and U components in the vicinity of the peculiar point and the randomly distributed noise. From Eq.(|29|) we 
can estimate the variance of the random function /: 
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(30) 



which depends on the module P{x,y). From Eq.(^) we can define more precisely the limits of applicability of the 
approximation Eg. (|29|) : > af^ and therefore 

an < ^P. (31) 

This relation shows that influence of the noise on the contour lines around peculiar point of polarization field becomes 
stronger when we come closer to the point P — 0. 

Finally at the very vicinity of the point P = we come to the region where > ■^P and for this region we have 
the approximation 

P^=nl + nl. (32) 



If uq and nu are independent random Gaussian functions then P = Uq + is the Rayleigh random process with 
pdf $([/): 

$([/) = -^expf-i^l. (33) 



As we can see from Eq.(B3|) the probability to find the point with zero level of polarization inside the area > 

\P [x, y) is proportional to P / (4tT^) and goes to zero if P ^ 0. This means that pixel noise or (and) polarized point 
sources destroy the structure of the field around nonpolarized points of the primordial signal and in principle remove 
this peculiarities from the map. Thus the regions around the peculiar points of polarization are very sensitive to the 
pixel and point sources noises. The structure of the polarization filed around such points gives us the information 
about the variance and statistical properties of the noises and can be used for filtration of such noises. We discuss 
this problem in details in [11] for filtration of the polarization maps. 



V. CONCLUSIONS 



New theoretical methods have been developed recently to provide the analysis and exploiting the new and forcoming 
CMB data sets. The BOOMERANG and MAXIMA-1 missions have discovered a lot of "hot" spots of the CMB 
radiation on the sky. These peaks in AT-distribution on the sky play a role of "standard" sources of the CMB 
anisotropy signal for future balloon and satellite experiments with more sensitive radiometers and higher angular 
resolutions. The new theoretical methods are developing now to provide the analysis and predictions the new and 
forthcoming CMB data sets. 
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